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ABSTRACT 

Aims. We explore the use of the color-log(n) (where n is the global Sersic index) plane as a tool for subdividing the galaxy population 
in a physically-motivated manner out to redshift unity. We thereby aim to quantify surface brightness evolution by color-log(n) type, 
accounting separately for the specific selection and measurement biases against each. 

Methods. We construct (u-r) color-log(n) diagrams for distant galaxies in the Hubble Ultra Deep Field (UDF) within a series of 
volume-limited samples to z = 1.5. The color-log(n) distributions of these high redshift galaxies are compared against that measured 
for nearby galaxies in the Millennium Galaxy Catalogue (MGC), as well as to the results of visual morphological classification. Based 
on this analysis we divide our sample into three color-structure classes. Namely, 'red, compact', 'blue, diffuse' and 'blue, compact'. 
Luminosity-size diagrams are constructed for members of the two largest classes ('red, compact' and 'blue, diffuse'), both in the UDF 
and the MGC. Artificial galaxy simulations (for systems with exponential and de Vaucouleurs profile shapes alternately) are used to 
identify 'bias-free' regions of the luminosity-size plane in which galaxies are detected with high completeness, and their fluxes and 
sizes recovered with minimal surface brightness-dependent biases. Galaxy evolution is quantified via comparison of the low and high 
redshift luminosity-size relations within these 'bias-free' regions. 

Results. We confirm the correlation between color-log(n) plane position and visual morphological type observed locally and in other 
high redshift studies in the color and/or structure domain. The combined effects of observational uncertainties, the morphological K- 
correction and cosmic variance preclude a robust statistical comparison of the shape of the MGC and UDF color-log(n) distributions. 
However, in the interval 0.75<z<1.0 where the UDF /-band samples close to rest-frame S-band light (i.e., the morphological K- 
correction between our samples is negligible) we are able to present tentative evidence of bimodality, albiet for a very small sample 
size (17 galaxies). Our unique approach to quantifying selection and measurement biases in the luminosity-size plane highlights the 
need to consider errors in the recovery of both magnitudes and sizes, and their dependence on profile shape. Motivated by these results 
we divide our sample into the three color-structure classes mentioned above and quantify luminosity-size evolution by galaxy type. 
Specifically, we detect decreases in S-band, surface brightness of 1 .57 ± 0.22 mag arcsec"^ and 1 .65 + 0.22 mag arcsec"" for our 'blue, 
diffuse' and 'red, compact' classes respectively between redshift unity and the present day. 

Key words, galaxies: evolution - galaxies: formation - galaxies: high-redshift - galaxies: photometry 



^ 1. Introduction 



Recent results from wide field imaging surveys, such as the 



Sloan Digital Sky Survey (SDSS; Adelman-Mc Carthv et al 



20061) and the Millennium Galaxy Catalogue (MGC dUske et al 



2003h . have highlighted the bimodal nature of the local galaxy 
population. In particular, bright galaxies at z ~ natu rally divide 
into two distinct types in the color-m agnitude ( Strateva et al.' 
12001 ). log(»)-m agnitude ([Driver et al.1 12006: BaUetal. 200^ 
and color-log(«) ([Driver et al. ||2006|) planelTOne type displays 
red colors consistent with old stellar populations (and/or very 
high dust contents) and the compact light profiles typical of 
pressure-supported systems. The other type has bluer colors 
indicative of more recent star-formation and younger stellar 
populations, and t he diffuse l i ght pr ofiles of rotationally- 
supported systems. [Driver et"aLl ( 20061) demonstrate that the 
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' Here n refers to the global 'Sersic index', i.e., the inverse-exponent 
in the best-fit Sersic model to the galaxy light distribution. The Sersic 
model is "most commonly expressed as an intensity profile, such that 



m 



,exp 



1 >, where is the intensity at the ef- 



fective radius Re that en closes half of the total light from the model" 
dGraham & Drivedl2005h and b„ is a normalization constant. 



red, compact peak generally corresponds to morphologically 
early-type systems (E/SOs), while the blue, diffuse peak is 
dominated by late-type systems (Sd/Irr), with intermediate 
spirals (Sa-Sc) spanning both. This bimodality suggests that 
there are two distinct pathways of galaxy evolution and for- 
mation. It is incumbent upon theories of galaxy formation to 
explain the origin and nature of these two pathways, and the 
evolution of galaxies within and between them. In the orthodox 
hierarchical clustering (HC) scenario the blue sequence (late- 
type population) is generally thought to form via the cooling 
of gas in rotating dark matter halos spun up by tidal torques 
in the early universe ( PeeblesllT969t [Fall & Efstathiou[[l98a) . 
The red sequence (early-type population) i s then built up via 
mergers between blue sequence members ( [Toomre & Toomr3 
1972; Toomre 1977 ) and dry mergers amongst the red sequence 
(Ivan Dokkumet all [1999[: [Khochfar & BurkertI [2003h . Fading 
of outer discs in spirals to reveal a central bulge, or conversely 
enhancement of the bulge relative to the disc, may also be 
an important evolutionary pathway ([Christlein & Zabludofj 
20i)4 and references therein). Secular evolution, such as 
the formation of bars and b ulges via disc instabilities (see 
[Kormendv & KennicuttI [2004 1 for a rece nt review) and perhap s 
the formation of outer disc truncations ([Debattista et al.[|2006h . 



2 



E. Cameron and S. P. Driver: Galaxy evolution by color-log(n) type since z ~ 1 in the UDF 



will likely play a significant role in isolated systems on long 
timescales. 

Deep imaging surveys, such as the Hubble Ultra Deep Field 
(UDF), allow us to constrain quantitative predictions from 
galaxy formation theories and gauge the importance of specific 
evolutionary mechani sms at different epochs. For instance, 
IMcIntosh et al] (l2005h compare high redshift galaxies in the 
GEMS survey to those at low redshift in the SDSS and find 
that the luminosity-size and mass-size relations are consistent 
with passive ageing of old stellar populations since z ~ 1. 
By a further comparison to theoretical models they limit the 
role of dissipationless merging to <1 major merger over this 
timescale, favoring instead the role of disc fading in order 
to build up the red sequence number density. Alternatively, 
[Blanton (2006) examines evolution of the luminosity functions 
of the red and blue sequences from z ~ 1 in DEEP2 to z ~ in 
the SDSS and finds negligible change in the number density on 
the blue sequence (to within 10%). He argues that this supports 
the need for dry mergers to evolve the red sequence number 
density given the lack of luminous blue progenitors. The crucial 
elements of this genre of study are the system adopted for 
assigning galaxy type and the methodology employed to assess 
selection and measurement biases. Recent work has favored 
physically-motivated, objective and quantitative classification 
systems for the former For instance, one can make use of 
the correlation betw een morphological type a nd Sersic index 
( Tru)illo et al.ll2006 l) or color and Sersic index (IMcIntosh et al.l 
l2005h to separate early-type and late-type galaxies. In order 
to maximize the reliability of the division between physical 
types one might attempt to identify each peak of the bimodality 
and place the division across the corresponding saddle point 
at each epoch. This approach was successfully applied by 
iBell et al.l d"2004 and [Blanton (2006) studying the high redshift 
color-magnitude bimodality, but has yet to be applied to the 
joint color-log(n) bimodality. 

Understanding the nature of the selection effects that operate 
on any given imaging survey is crucial to the robust interpre- 
tation of the data. These include biases against the detection 
and photometry of compact, extended, faint and/or low surface 
brig htness systems . This is well-documente d in the literature 
(e.g. iDis nev 1976; Disnev & Phillipps 1983; Impev & Bothuri 
119971; [Driver 1999; Cross & Driver 2002; Driver etal. 2005). 
The low surface brightness limit is of vital importance to studies 
that attempt to recover the properties of apparently faint and/or 
dim galaxies, whether at the faint end of the luminosity function 
at low redshift or for intrinsically bright galaxies at high red- 
shift. As Dalcanton (1998) and Cross & Driver (2002) point out, 
low surface brightness systems suffer from both a Malmquist 
incompleteness bias and from an under-estimation of their 
total fluxes — and the dual impact of these sources of error can 
conspire to mimi c the effect of gala xy evolution with redshift . 
For instance, Simard et alj (1 19991) and lRavindranath et all (|2004 
argue that the surface brightness evolution reported in the first 
generation of Hubble Space Telescope (HST) imaging survey 
studie s (e.g. Schade et al. 1995; Liflv et al. 1998; Roc he et all 
[1991 was almost entirely the result of sel ection bias. More re- 
cent, deep, HST, imagin g studies (such as IMcIntosh et al.|[2005[ 
and iTruiillo et alj l2006l) have adopted improved approaches 
to handling selection bias (and do, in fact, find evidence for 
a significant degree of surface brightness evolution to z ~ 1). 
However, best practice methods are not applied universally 
and there is as yet no clear consensus on the best way to 



evaluate surface brightness bias. In ICameron & Drivel (|2007[) 
we introduced a relatively conservative approach for identifying 
bias-free regions in the luminosity-size plane. We demonstrated 
that even in the exceptionally deep UDF imaging the systematic 
under-estimation of total flux by Kron-like aperture magnitudes 
at low surface brightness could be causing bright, yet diffuse, 
systems to be recovered as faint, compact objects. However, in 
that study we limited our investigation to the case of exponen- 
tial, disc-like profiles, but for mor e concentrated galaxies the 
biases are known to be much worse dGraham & Driveill2005h . It 
is, therefore, important to extend this analysis to 'elliptical-like', 
de Vaucouleurs profile shapes in order to properly quantify 
evolution in the 'red, compact' population. 

In this paper we construct the color-log(n) distribution of 
galaxies out to z ~ 1.5 in the Hubble Ultra Deep Field (UDF) 
and compare it against that of the local galaxy population in 
the Millennium Galaxy Catalogue (MGC). We examine the 
use of this parameter space for morphological classification 
and search for evidence of bimodality. Dividing our samples 
into 'blue, compact', 'blue, diffuse' and 'red, compact' classes, 
we then investigate evolution in luminosity-size relations of 
the latter two (predominant) classes. Careful attention is paid 
to the relevant selection effects and uncertainties, which we 
evaluate via artificial galaxy simulations. The outline of the 
paper is as follows: in Section 2 we describe the construction 
of our dataset; in Section 3 we quantify selection effects and 
measurement biases; in Section 4 we present our results and 
discuss the relevant uncertainties; in Section 5 we compare our 
findings to other recent studies and discuss their interpretation 
in light of the hierarchical formation paradigm; and in Section 
6 we summarise our conclusions. A cosmological model with 
Q.M = 0.3, = 0.7 and Hq = 100 km s"' Mpc"' is used 
throughout. These specific values of the cosmological parame- 
ters were adopted for ease of comparison between the present 
UDF work and the slightly older MGC results. Unless otherwise 
stated, all magnitudes are given in the AB system. 



2. Data 

The Hubble Ultra Deep Field d Beckwifli et al.l l2006l) consists 
of an 11 arcmin^ patch of sky centred on RA = 03° 32' 39.0", 
Dec = -27° 47' 29.1" (J2000) in the region of the Fornax 
Constellation. The publicly released ACS/WFC Combined 
Images (Version 1.0) span the optical wavelength range 3700 
to 10,000A in four wide-band filters: F435W (B), F606W 
(V), F775W (0 and F850LP (z). The F775W /-band image 
has the longest total exposure time of 347,110 s (144 orbits). 
Each single exposure was half an orbit in duration with the 
pointing cycled through a four part dither pattern. Each image 
has been processed through the standard HST data pipeline 
and drizzled to a pixel scale of 0. 03" pixel '. The PSF FWHM 
is ~0."081 in the i-band image dBeckwifli et al.ll2006l) . which 
is close to the theoretical diffraction limit of the HST at these 
wavelengths. Supplementary observations in the near-IR have 
been made wifli NICMOS in the FllOW (J) and F160W {H) 
filters dThompson et alj[2005h . although only for a reduced area 
of ~5.76 arcmin^ within the UDF region. These images are 
drizzled to a coarser pixel scale (0."09 pixel ') appropriate to 
the larger PSF of the NIC 3 camera (~0."36 FWHM in flie 
7-band). 
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2.1. The iUDF-BRIGHT Sample, Kron Magnitudes and Half 
Light Radii 

In ICameron & Driverf (l2007 h we presented iUDF-BRIGHT— a 
catalogue of 2497 sources brighter than 28th magnitude in the 
ACS i-band image. The se were identified using SExtractor 
jBertin & ArnoutsI Il996h with a constant surface brightness 
threshold of 27.395 mag arcsec"^ and were matched, where 
possible, with counterparts in the Version 1 .0 data release (BOfl) 
catalogue to compare object magnitudes and deblending. Our 
isophotal magnitudes were found to be in good agreement with 
those in B04 with the largest mean difference just 0.06 + 0.07 
mag in the worst 27-28th magnitude bin. A total of 125 objects 
were found to have conflicting centroid positions, so these were 
visually inspected. The B04 deblending was preferred in 60 
instances, our deblending preferred in 56 instances and there 
were 9 spurious detections. Elliptical aperture, major axis, half 
light radii were then computed for the entire iUDF-BRIGHT 
sample using Kron magnitude fluxes and a crude PSF correction. 



2.2. Colors, Photometric Redshifts and Sersic Indices 

ICoe et all (l2006l) present magnitudes, photometric redshifts and 
Sersic indices for a large sample of objects in the ACS and 
NICMOS UDF images. Their catalogue (C06) includes both 
the B04 / and z-band sources and the T05 J+H-hand sources, 
as well as their own new B+V+i+z and J +H detections. The 
magnitudes provided by ICoe et all ( l2006l) for each filter are 
well suited for computing galaxy colors as they are aperture- 
matched, PSF-corrected, and corrected for the difference 
between isophota l and Kron fl uxes. Photometric redshifts are 
estimated by Coe et al.l (l2006l) using the Bayesi an BPZ code 
(lBenftedl200ar They adopt the SED library of iBenitez et alj 
(12004') with the addition of two new templates representing 
young, blue, star-forming galaxies. Their photometric redshifts 
agree well with those determined spectroscopically for 41 bright 
galaxies in the UDF with an RMS of Az - 0.04(1 + Zspec) 
up on exclusion of 4 outlier s. However, as we demonstrated 
in ICameron & Drivel (12007 ^ by comparison to an alternative 
UDF photometric redshift catalogue (M05, supplied by B. 
Mobasher, priv. comm.), the redshift estimates derived using 
this technique are highly sensitive to the choice of SED library 
and input magnitude type, particularly for faint galaxies. This 
contributed a large source of uncertainty to our previous results, 
although we restrict our investigation to bright galaxies with 
Mb < -18 mag and 0.25 < z < 1.5 for which the redshifts 
should be the most reliable. In fact, 80% of galaxies within these 
limits have their C06 and M05 photometric redshift es t imates 
in agreement within +0.1(1 -i- Zav)- Finally, ICoe et alj (l2006l) 
provide Sersic indices for all objects in the B04 catalogue. 
These are global Sersic indices computed via 2D model fits 
to each galaxy's /-band image using the galfit software package. 

As mentioned above, in Cam eron & Drivel (l2007h we 
matched objects in the iUDF-BRlGHT sample to those in 
B04 to check o ur source extraction and photometry. Since 
ICoe etalJ ([20061) preserve B04's ID numbering system, it was a 
simple process to link the C06 morphological and photometric 
catalogues to our galaxy catalogue. None of the 56 galaxies for 
which we preferred an alternative deblending to that of B04 

2 In|c ameron & Drive j ( l2007|) we called this 'the on-line catalogue', 
but here we adopt t he syntax of ICoe et al. ( 2006) for simplicity. Hence, 
we also now refer to lThompson et alj l2005,) 's J+H catalogue as 'T05' . 



have redshifts in C06 that place them within the limits of this 
study. Hence, their lack of a proper match does not affect our 
work. 



2.3. K-Corrections, Absolute Magnitudes and Rest-Frame 
Colors 

Individual galaxy K-corrections are needed to compute rest- 
frame absolute magnitud es and colors. Here we follow a similar 
procedure to ICameron & Driver (2007), fitting redshifted 
spectral templates from the library of Poggianti d 19971) to 
each galaxy's observed optical and near-IR (where available) 
broadband photometry. However, this time we use the aperture- 
matched, PSF-corrected fluxes and BPZ redshifts from the C06 
catalogue. The magnitudes used in our luminosity-size relations 
are all derived by K-correcting from the observed /-band (to 
rest-frame MGC B-band) because we have performed detailed 
modelling of the detection completeness and photometric 
reliability of galaxies in the ACS /-band image. However, when 
deriving rest-frame (u - r) colors we utilize bandpass shifting, 
i.e., K-correcting from the observed filters sampling nearest 
to rest-frame u and r-band light at each redshift (either side 
of the 4000A break). Hence, we correct to (m - r)iest from 
{B - /)obs at 0.25 < z < 0.5, (V - z)obs at 0.5 < z < 1.0 and 
(/ - 7)obs (if available, otherwise (V - z)obs) at 1.0 < z < 1.5. 
Absolute magnitudes and colors are then derived using these 
K-corrections and a cosmological model with Qm - 0.3, 
Qa = 0.7 and Hq - 100 km s"' Mpc These specific values of 
the cosmological parameters were adopted in order to facilitate 
comparison to the earlier, published MGC data. 



3. Selection Effects 

Here we quantify the selection effects and systematic measure- 
ment biases for galaxies in our sample in the luminosity-size 
plane. Our approach is more conser vative than o ther recent 
studies (e.g. Bouwens et al.ll2004t ISargent et al.ll2007l) because 
we use artificial galaxy simulations to test not only object 
completeness, but size and magnitude recovery as well. We 
then select only the region of the luminosity-size plane at each 
redshift where the detection completeness of our simulated 
galaxies is over 90% and at least 75% of these have measured 
fluxes and half light radii within 10% of their input (i.e., 
intrinsic) values. We previously applied this methodology 
successfully to our study of the z ~ 0.75 UDF ga laxy population 
(not subdivided by galaxy type) in ICameron & Driver (2007'), 
but only using exponential {n - 1) profile shapes (on the 
assumption that our sample was dominated by disc-dominated 
galaxies at this redshift). Here we extend this work to cover 
compact, bulge-dominated galaxies by additionally performing 
simulations using de Vaucouleurs (« = 4) profile shapes, and 
contrasting the results against those from our exponential profile 
simulations. 

For our UDF galaxy photometry we use SExtra ctor's 
Kron magnitudes and half light radii. Kron magnitudes dKronl 
Il980l) employ a luminosity-weighted average to define the 
so-called Kron radius {R\), and the corresponding magnitude 
is computed by summing the flux within an aperture of radius 
2.5R\. Theoretically, the Kron magnitude should contain 96.0% 
of the flux in an exponential profile and 90.4% of the flux in 
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a de Vaucouleurs profile (integrated to infinity). However, in 
practice these theoretical values are rarely obtained because of 
observational constraints. As Graham & Driver (2005) point 
out 'what is important is not the sky-level or isophotal-level 
one reaches, but the number of effective radii that have been 
sampled' . It is this source of error that we are effectively trying 
to gauge via our artificial galaxy simulations. Of course, when 
we refer to our ability to accurately recover the total or true flux 
of a de Vaucouleurs type profile we mean the total Kron flux 
and are allowing for the 9.6% loss against the integration of 
the profile to infinity. As Kron magnitudes are used in both the 
UDF and MGC photometry this gap between the total integrated 
and total Kron magnitudes will not introduce a source of error 
into comparisons of their relative magnitude, size and surface 
brightness distributions. It is also important to note that errors 
in the measured Kron magnitude will flow on to cause errors 
in the measured half light radius, with an under-estimation of 
flux leading to an under-estimation of size. The extent of this 
under-estimation again varies with profile shape. Hence, we test 
our ability to accurately recover both magnitudes and sizes here. 



3.1. Artificial Gaiaxy Simuiations and Apparent Selection 
Limits 

Our procedure for performing these selectio n bias tests is 
outlined in detail in ICameron & Drivej (l2007l) . Essentially, it 
involves dividing the observed luminosity-size plane into a grid 
and at each grid point generating 100 artificial galaxies with 
random inclinations and axial ratios using the IRAF artdata 
package. These are inserted into the UDF i-band image, and 
SExTRACTOR used to search for them and to measure Kron mag- 
nitudes. Half light radii are then computed using the elliptical 
aperture method. The resulting detection positions, magnitudes 
and sizes are compared to those in the input object list. 

Figure [1] displays the results of our simulations for galaxies 
with de Vaucouleurs (n = 4) profile shapes in a set of three 
plots. The first plot indicates our detection completeness as a 
function of apparent magnitude and effective radius (i.e., the 
number of artificial galaxies detected of the 100 inserted at 
each grid point). The second plot reveals our 'reliability' in this 
parameter space, which we define to be the number of detected 
galaxies having measured magnitudes and effective radii within 
10 per cent of their input values. And finally, we produce an 
'error vector' diagram showing the (3cr-clipped) mean size and 
direction of the difference between input and recovered values 
for the detected objects originating in each bin. Together, the 
reliability and error vector diagrams allow us to identify any 
regions of the observable parameter space contaminated by 
galaxies with measured magnitudes and half light radii that 
poorly reflect their true, intrinsic properties. Equivalent plots for 
ou r exponential (n - 1) pro file shapes were already presented 
in ICameron & Drivel (12007) . although here we have tightened 
the tolerance on flux and magnitude recovery from 25% to 10%. 
Thus, we recompiled our previous n - 1 profile results with 
this new tolerance so as to be able to indicate the corresponding 
75% completeness and reliability limits for this profile shape for 
comparison in Fig.[T] 

It is clear from our completeness diagram that de 
Vaucouleurs (n - 4) profile shape galaxies are easily detectable 
over the entire area of the apparent magnitude-half light radius 
plane spanned by the observed iUDF-BRlGHT distribution. 
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Fig. 1. Completeness, 'reliability' and 'error vector' diagrams 
for de Vaucouleurs (n = 4) profile objects in the ACS UDF i- 
band image computed from artificial galaxy simulations as de- 
scribed in Section 3.1. The greyscale squares and the arrows in- 
dicate the simulation results and the black dots show the dis- 
tribution of the real observed data. The 75% completeness and 
reliability limits that define our bias-free region are marked for 
clarity. The equivalent limits for exponential (n = 1) type galax- 
ies at low surface brightness are also indicated for comparison. 
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Their highly centrally-concentrated light distributions ensure 
that even for low mean surface brightness systems there is still 
a compact, high surface brightness core that emerges above the 
background sky. Hence, they are more readily detectable than 
exponential profile galaxies which we fail to identify in the 
UDF /-band image beyond ~28.0 mag arcsec"^. However, when 
we investigate how accurately we can measure the intrinsic 
magnitudes and sizes of these galaxies the situation is reversed. 
It is actually extremely difficult to recover the intrinsic flux of 
a de Vaucouleurs profile system as a significant fraction of the 
light is contained in its low surface brightness outer regions (the 
profile 'wings'). Whilst we can recover intrinsic magnitudes 
and sizes of exponential profile galaxies to within 10% accuracy 
down to {n)g~21 .0 mag arcsec^^, here we find we can only 
recover those of de Vaucouleurs profile galaxies down to ~26.0 
mag arcsec"^. So one might wonder, ;/ there exists a significant 
population of high Sersic index systems intrinsically occupying 
the region of the apparent magnitude-size plane in between our 
completeness and reliability limits, where do they turn up in 
the observed distribution? For this we inspect the 'error vector' 
diagram, which indicates they are erroneously recovered as 
faint, compact objects. This illustrates the importance of using 
a surface brightness selection limit based on measurement 
reUability as well as detection completeness. 

The reliability plot also reveals constraints on our ability 
to accurately observe very large, very small and very faint 
galaxies. The properties of galaxies larger than ~4" in the UDF 
are difficult to measure as they overlap with numerous other 
sources, causing problems with deblending. Galaxies smaller 
than ~0".06 pose a challenge for accurate measurements of the 
half light radius as they are a similar size to the PSF (~0".081) 
and pixel scale (0".03). Also, we cannot measure reliable 
magnitudes for objects fainter than 28.0 mag in the /-band 
image. This is because, even though 28.0 mag is well above the 
lOo" point source detection limit of 29.2 mag, the signal-to-noise 
ratio is low outside the cores of faint, extended sources. Due to 
the decision to locate the UDF away from known, nearby, bright 
objects there is also an effective selection limit on the brightest 
galaxies in the sample. One can estimate this limit from the 
brightest apparent (galaxy) magnitude in the catalogue, which is 
18.26 mag in the /-band. 

By considering all the limits described above we define 
an observational window in the apparent magnitude-size plane 
inside which all galaxies are detected with high completeness 
and structural parameters are reliably recovered. We shall refer 
to this as the 'bias-free' region. This five-sided, bias-free region 
is indicated in Fig.[T| 



3.2. Absolute Selection Limits 

The selection boundaries derived for the apparent luminosity- 
size plane are easily extended to the absol ute reg i me fo r 
a volume-limited sample via the method of iDriveii (Il999h . 
The relevant eq ua tions are presented and explained in 
ICameron & Driven (l2007b . and include the effects of both 
the (1 + z)'^ cosmological surface brightness dimming and the 
K-correction. In this study we use a series of 5 redshift bins 
for our volume-limited samples spanning the redshift range 
z - 0.25 to 1.5. Our absolute selection limits for galaxies in 
each redshift bin are marked on the luminosity-size diagrams of 
Fig.|6] Note that we apply the « = 1 limits to our 'blue, diffuse' 



class and the n - 4 limits to our 'red, compact' class. 



4. Results 

4.1. Color-Log(n) Diagrams 

In Fig. |2] we present the iUDF-BRIGHT color-log(n) distribu- 
tion of bright galaxies (Mb < -18 mag) in five volume-limited 
samples from z = 0.25 to 1.5. We also plot the color-log(n) 
distribution of the 7877 bright galaxies {Mb < -18 mag) at 
0.0 < z < 0.25 in the Millennium Galaxy Catalogue (MGC) 
(iLiske et al.l l2003t [Driver et al.l 120051) fo r comparison as a low 



redshift benchmark. Driver et al. (120061) demonstrate that the 



MGC local galaxy sample shows clear bimodalities in both 
rest-frame (m - r) color and global Sersic index, and that 
there exist two extremely well separated populations in the 
color-log(n) plane. This bimodality in the benchmark MGC 
distribution is highlighted in Fig. |2]by volume-density contours 
overlaid against the raw MGC distribution in the first panel and 
against the iUDF-BRIGHT high redshift samples in the subse- 
quent panels. The saddle point between the two peaks offers 
an objective, physically-motivated subdivision of the galaxy 
population. We hoped to detect the equivalent bimodality in our 
high redshift bins, thereby ensuring self-consistent classification 
over all epochs. However, even with the MGC contours guiding 
the eye, no clear bimodality is evident in the iUDF-BRIGHT 
samples, except perhaps at 0.75 < z < 1 .0. This is presumably 
due to the combined effects of cosmic variance, luminosity 
evolution, measurement uncertainties and the morphological 
K-correction as we discuss below. The final divisions we use 
to separate our sample into 'blue, diffuse', 'red, compact' and 
'blue, compact' classes are indicated by thick black lines in 
Fig.|2] The choice of these divisions is explained in Section |4.2.2| 



4.2. Sources of Uncertainty and Bias 

It is important to bear in mind that the result s presented here 
are highly susceptible to cosmic variance (see ISomerville et"aLl 
|2004). Given the UDF field of view of 11 arcmin^ we note 
that the volumes surveyed within the various redshift intervals 
shown in Fig. 2 are: 608 Mpc^ (0.25 < z < 0.50); 1261 Mpc^ 
(0.50 < z < 0.75); 1860 Mpc^^ (0.75 < z < 1.00); 2332 Mpc^ 
(1.00 < z < 0.25); and 2674 Mpc^ (1.25 < z < 1.50). Within 
suc h small volumes cosmi c variance can be severe (up to 200%, 
c.f. lSomerville et aLll2004 Fig. 2) and any trends could therefore 
be due to both evolutionary and environmental effects. Based 
on random sampling of the local MGC survey one typically 
requires a volume of order 100,000 Mpc^ to reduce cosmic 
variance to ~10%, i.e., ~50 UDFs. 

One of the weaknesses of the color-log(n) bimodality as 
a tool for identifying galaxy types is its sensitivity to the 
magnitude Umit imposed. The original Mb < -18 mag cut was 
chosen as it gives roughly equal numbers of red and blue peak 
galaxies in the local MGC sample. With a brighter limit the red 
galaxies dominate the sample and with a fainter l imit the blue 



galaxies dom inate. Nume r ous studies (recen t ly iBarden et al 



2005; Dahlen et all 120051; lllbert et 



es (recen tly 
ID 120061: H 



Mcintosh et al 



120051: iZucca et al.ll2006l) have found evidence for a decrease in 
the characteristic luminosities of different galaxy types from 
z ~ 1.5 to the present, attributed primarily to the ageing of 
their stellar populations. Hence, our fixed absolute magnitude 
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Fig. 2. Color-log(n) diagrams for the MGC 0.0 < z < 0.25 sample and t he iUDF-BRJGHT volume-hmited samples to z = 1.5. 
The MGC local galaxy distribution is clearly bimodal (as demonstrated bv lDriver et al.ll2006l) and its number density contours are 
overlaid on the high redshift UDF panels for comparison as a local benchmark. The UDF distributions are not evidently bimodal, 
except perhaps at 0.75 < z < 1.0, for reasons discussed in Section l4!2l including cosmic variance and measurement errors. Eyeball 
morphological classifications for iUDF-BRIGHT galaxies at z < 1.0 are indicated with symbols corresponding to solid square 
= E/SO, solid triangle - Sa-Sc and open hexagon = Sd/Irr/Merger The solid lines indicate our division of the sample into three 
categories, namely 'blue diffuse', 'red compact' and 'blue compact'. Objects in all panels are subject to a constant faint magnitude 
limit of Mb < -18 mag. We also explore the effect of brightening this limit to account for the passive ageing of galaxy stellar 
popuations (see Section l4l2l i. Objects fainter than the evolving magnitude Umit are colored grey instead of black. 



limit of Mb < -18 mag may include some galaxies at high 
redshift that would become too faint over time to make the 
cut at low redshift. We attempt to account for this effect by 
correcting our high redshift hmits according to the anticipated 
luminosity evolution (as in iBarger et all 119981) . Specifically, 
we adopt corrections of AM = -Q.5z for our 'red, compact' 
class limit, AM = -0.6z for our 'blue, diffiise' class limit and 
AM - -1.75z fo r our 'blue, compact' clas s limit, ni o yitate d 
by the results of IZucca et al.1 ((2006) and filbert et al.1 (|2006|) . 
Any points excluded by these additional constraints are colored 
grey (instead of black) in Fig. |2] and Fig. |6] A large number of 
faint, blue galaxies are removed by this cut, but no significant 
increase in bimodality is evident. Interestingly, all but one of 
our genuine 'blue, compact' objects would likely be fainter than 
the Mb < -18.0 mag cut locally. It is likely these galaxies are 
the pre-cursors of modern day intermediate and dwarf galaxies 
caught at an epoch of intense star-formation. 

Measurement errors in individual Sersic indices and colors 
contribute only a small so urce of u ncertainty to the recovered 
UDF galaxy distributions. Coe et al.l 12006) estimate an error on 
the optimum model fit to each galaxy using galfit, which for 
our sample has a median value of An ~ 0.2. They also calculate 
uncertainties on their aperture magnitudes based on image 
noise, and the corresponding median values for our sample 
are AB = 0.01, AV = 0.06, Ai = 0.04 and Az = 0.05 mag. 
This implies median errors on our apparent (B - i) and {V - z) 
colors of 0.05 and 0.08 mag respectively. Additional errors of 
^0.04 mag are expected from our K-corrections given the SED 
template and band-pass shifting approach adopted. The total 
uncertainties on our rest-frame (u - r) colors are thus $0.09 mag. 

Finally, we note the problem of the morphological K- 
correction, i.e., the intrinsic variation in galaxy morphology 
with rest-frame wavelength. Coe et al. (2006)'s Sersic fits are 
all performed in the observed /-band, which at z ~ 0.75 samples 
close to rest-frame Z?-band (optical) light, but by z ~ 1 .5 samples 



closer to rest-frame J7-band (ultra-violet) light. The dominance 
of UV emission by bright, young stars is known to cause galaxy 
images to appear clu mpier and of later type than they do at 
optical wavelengths dKuchinski et al.l 120001: iTavlor-Mager et al] 
2007) and gal axy hght profile s to ap pear less concentrated. 
For instance, iLa Barbera et al.l (l2003l) find mean variations 
and standard deviations in Sersic index from optical to UV of 
<Alog(n)) = -0.09, 0- = 0.10 and <Alog(n)) = -0.02, cr = 0.19 
for samples of discs and spheroids respectively at z ~ 0.21. 
This effect may contribute towards obscuring any intrinsic 
bimodality in the color-log(n) distributions of UDF galaxies in 
the lowest and highest redshift bins of Fig. |2] 



4.2.1. Bimodality at 0.75 < z < 1.0 



Given the potential sources of uncertainty and bias that may 
distort our measured galaxy color-log(n) distributions at high 
redshift, it is perhaps no surprise that most of the UDF volume- 
limited samples fail to reveal strong evidence of bimodality. 
However, in the 0.75 < z < 1.0 bin where the morphological 
K-correction is negligible there does appear to exist two distinct 
(well separated) populations — one displaying the blue colors of 
star-forming galaxies and the difiiise light profiles (low Sersic 
indices) of disc-dominated systems, and the other displaying 
the red colors of old stellar populations and the compact light 
profiles (high Sersic indices) of bulge-dominated systems. 
However, as mentioned earlier, the small volume sampled by 
the UDF over this redshift interval leaves this result highly 
susceptible to cosmic variance. Further investigation of this 
issue will thus be warranted once a suitably large, deep imaging 
dataset becomes publicly available. 
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4.2.2. Sample Division and IVIorpliological Correlation 



4.3. Quantifying Color Evoiution 



Our main aim in constructing the UDF high redshift color- 
log(«) diagrams was to evaluate their suitability as a tool for 
subdividing the galaxy population. Unfortunately, the lack of 
clear evidence of bimodality in these high redshift distributions 
means that our classification approach cannot be defined 
self-consistently from the galaxy populations at each epoch. 
Hence, we adopt the saddle point of the local bimodality, 
(u - r)i-est = 2.325 - 2.074 log(«), to separate the two key 
galaxy types, i.e., our 'blue, diffuse' and 'red, compact' classes. 
Visual inspection of the data reveals that this is still an appro- 
priate division of the two distinct populations apparent in the 
0.75 < z < 1 .0 bin. We note also note that the steepness of this 
line in color-log(n) space reduces the sensitivity of our cut to 
the observed evolution in galaxy colors. Postage stamp images 
showing examples of each of these galaxy classes are shown in 
Fig. [3] and Fig.|4] In addition we classify galaxies in the region 
described by log(n) > 0.398 and (u - r)rest < 1.5 as 'blue, 
compact' systems. These systems are common in our sample 
at high redshift, but are unmatched by bright local counterparts 
in the MGC. Examples of this type of galaxy are shown in Fig. 
|5] They all exhibit a bright, compact nucleus embedded in an 
extended region of diffuse, irregular emission. Such galaxies 
have been observed in other surveys, predominantly at high 
redsh i ft (e.g.llm et al.| [2001; Menanteau et al. 2004; Cross et d] 
l2004t Illbert et al.1 l2 006t lEUis et all I2005t IDriver et all I2007Ik 
but also locally (e.g. [Driver et al. ri2005l). When the evolving 
luminosity limit was introduced, all but one of these galaxies 
were too faint for the cut. These galaxies are presumably not the 
direct pre-cursors of modern day giant systems, but are probably 
lower mass systems caught at an epoch of intense star-formation. 

Although we are attempting to implement an objective, 
quantitative method of classification, it is still informative to 
compare our cuts against the d istribution of visual morpho- 
logical types. [Driver et"aLl (l2006h reported a strong correlation 
between visual morphology and color-log(n) plot position. 
Namely, that early-types occupy the red, compact peak and 
late-type spirals and irregulars occupy the blue, diffuse peak 
with intermediate spirals bridging the gap. Hence, we inspected 
the ACS B, V, i and z-band postage stamp images of all UDF 
galaxies in our Mb < -18 mag sample with z < 1 and assigned 
them to one of three groups based on their Hubble type: E/SO, 
Sa-Sc, or Sd/Irr. Different symbols are used for each galaxy 
in our color-log(n) diagrams of Fig. |2] to indicate these types. 
Our high redshift galaxies echo the findings of [Driver et al] 
(1200 6) regarding the relation between color-log(n) position 
and visual morphology at low redshift. Furthermore, we see 
that our chosen division line between 'blue, compact' and 
'red, compact' types is indeed effective at separating E/SOs 
from Sa-Sc and Sd/Irrs. We do not attempt to assign eyeball 
morphologies to the 1.0 < z < 1.5 objects as these are more 
difficult to classify, being less well resolved in our images and 
dominated by irregular features. This is due both to intrinsic 
evolution in the star-formation rate and the fact that the / and 
z-bands are sampling rest-frame UV light (i.e., young stellar 
populations) at these redshifts. Even so, when we examine z > 1 
galaxies from each peak alongside their counterparts at z < 1 
there is a clear consistency in appearance, as seen in Fig. [3] and 
Fig- El Overall, this investigation of the visual morphologies of 
galaxies in our sample gives us confidence in our method of 
subdivision on the color-log(n) plane. 



A certain degree of color evolution is clearly evident in the 
color-log(n) distributions of Fig.|2] Specifically, the mean color 
of our 'blue, diffuse' galaxies becomes substantially redder from 
<(m - '■),-est) = 1.12 + 0.09 mag in the UDF at 0.75 < z < 1.0 and 
1.28 + 0.09 mag in the UDF at 0.5 < z < 0.75 to 1.59 + 0.01 
mag in the MGC at < z < 0.25. This represents a change 
of 0.62 + 0.10 mag from z - 1 to z = 0. The 'red, compact' 
galaxies also become redder on average, but not as strongly, 
evoluting from ((« - r)rest) - 2.21 + 0.09 mag in the UDF at 
0.75 < z < 1 .0 and 2.08 + 0.09 mag in the UDF at 0.5 < z < 0.75 
to 2.42 + 0.01 mag in the MGC at < z < 0.25. This represents 
a change of 0.37 + 0.10 ma g from z - 1 to z = 0. These results 
are in line with those of B lanto 3 (I2OO6) comparing objects at 
z ~ 1 in the DEEP2 survey to those found in the nearby universe 
in the SDSS. In that study blue sequence galaxies were found to 
be redder in " '(g - r) color (i.e., (g - r) color using a rest-frame 
of z = 0.1) by 0.3 mag by z ~ 0, while red sequence galaxies 
were redder by 0. 1 mag. In (g - r) color the redshifted 4000A 
break lies within the g-band whereas in (m - r)rest the two filters 
lie either side of it. This means that (m - r),est color is more 
sensitive to changes in the size of the 4000A break. Using 
our SED template fits and filter throughput functions we can 
compute ^'^{g - r) colors for our galaxies. In this way we find 
our (m - r)rest color changes are equivalent to ^'^{g - r) color 
changes of 0.31 ± 0.10 mag and 0.21 + 0.10 mag for our 'blue, 
diffuse' and 'red, compact' populations respectively. 



4.4. The Luminosity-Size Relation by Galaxy Type 

In Fig.|6]we present the B-band, iUDF-BRIGHT luminosity-size 
relation by galaxy type, i.e., for our 'red, compact' and 'blue, 
diffuse' classes, in our five volume-limited samples to z = 1.5. 
The relevant selection limits on completeness and measurement 
reliability for high and low n systems (as calculated in Section 
3) are marked with grey shading on each plot. Objects removed 
by our evolving magnitude limits are again colored grey instead 
of black. The equivalent MGC local galaxy relations were 
constructed by adopting the same color-log(n) criteria as for 
the UDF samples. An (inverse) observable volume-weighted 
correction was then applied to generate the MGC number 
density contours. These are overlaid on all the high redshift 
panels for comparison as a local benchmark. The 100 Mpc^ 
iso-volume contour, within which MGC selection biases are 
well constrained, is indicated with a thick black line. 

It is clear from these diagrams that the UDF survey's bias- 
free observational window covers a similar region of parameter 
space out to redshift unity as that which the MGC does locally. 
With the exception of the 0.75 < z < 1.0 'red, compact' bin, 
all our UDF luminosity-size distributions at z < 1 are well 
separated from the relevant surface brightness selection limits. 
Hence, we are confident that the sizes and magnitudes measured 
for these objects are free of surface brightness-dependent biases 
within a 10% tolerance. At these redshift intervals the 'blue, 
diffuse' UDF galaxies appear to show evidence of evolution 
towards higher Z?-band surface brightnesses relative to the z ~ 
MGC population. Quantitatively, we measure mean (ju)^ of 
21.3 + 0.2 and 21.1 ± 0.2 mag arcsec"^ for galaxies brighter 
than the evolving magnitude limit in the UDF 0.5 < z < 0.75 
and 0.75 < z < 1.0 bins respectively. This is compared to 
22.20 + 0.01 mag arcsec"^ in the MGC. (The error margins are 
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Fig. 3. Galaxy morphology in our 'blue, diffuse' class. Each row contains postage stamp images showing members of this class cut 
from the ACS/WFC /-band UDF image. These are clearly disc-dominated systems, predominantly late-type spirals and irregulars. 
The rows are ordered in increasing redshift intervals from 0.25 < z < 0.50 to 1.25 < z < 1.5 to illustrate the changes in galaxy 
appearance due to bandpass shifting, decreasing angular sizes and morphological evolution. 



derived from the sum in quadrature of the uncertainties on the 
apparent magnitudes, sizes and redshifts for an individual UDF 
galaxy). This is a 1.57 + 0.22 mag arcsec"^ increase from z = 
to z = 1 (based on the line of best fit between the one MGC and 
two high redshift UDF datapoints) and is in good agreement 
with recent results from other high redshift studies. For instance, 
iBarden et all (12005 ) find evolution of 1.43 + 0.07 mag arcsec ' 
to z = 1 in rest-frame Z?-band for an My < -20 mag disc 
galaxy sample in the GEMS survey data. The slightly smaller 
amount of evolution they find could be due to the brighter limit 
of their sample as there is some evidence for a luminosity (or 
size) dependence to the surface brightness evolution in our 
data. Least squares fits to the luminosity-size relations in the 
UDF 0.5 < z < 0.75 and 0.75 < z < 1.0 bins return gradients 
of -0.21 + 0.06 and -0.24 + 0.06 log(kpc)/mag respectively, 
compared to -0.177 + 0.004 log(kpc)/mag for the MGC local 
relation. Other st udies have yielded similar results. For instance, 
iMelboume et al.l ([2007 ) report Z?-band luminosity evolution of 
1.53+q [q mag to z = 1 for large (Re > 5 kpc) blue galaxies in 
GOODS-N and 1.65+jJ 'j'^ mag for intermediate (3 < Re < 5 kpc) 
ones. Beyond z ~ 1 our surface brightness-dependent limits on 
UDF object detection and parameter reliability impinge on the 
observed distribution and the valid region of comparison to the 
MGC data. Thus, although the mean (/i),, continues to increase 
to ~20.7 mag arcsec"^ by our 1.0 < z < 1.25 bin, we cannot 
distinguish between real evolution and the impact of selection 



effects at these redshifts. 

The high redshift luminosity-size relations for our 'red, 
compact' sample are somewhat more difficult to interpret. Even 
though the selection limits are significantly tighter on these 
higher Sersic index systems (as shown in Section 3.1), our 
UDF 0.5 < z < 0.75 bin still samples an equivalent region 
of parameter space to the MGC survey, and our galaxies are 
detected and measured free of bias. One can see from Fig.|6]that 
the luminosity-size relation of the local MGC 'red, compact' 
population has a higher surface brightness and steeper slope 
than the 'blue, diffuse' one, which partially compensates for 
the tighter selection limits. This is also the case in our UDF 
0.5 < z < 0.75 bin where the mean 'red, compact' is 
20.1 + 0.2 mag arcsec"^, i.e., 1.2 mag arcsec"^ brighter than 
the 'blue, diffuse' galaxies at this redshift. It is also an increase 
in surface brightness compared to the MGC local sample of 
'red, compact' systems, which has a mean (yu)^ of 21.65 + 0.01 
mag arcsec"^. Extrapolating this trend with a linear fit gives an 
evolution of ~3 mag arcsec"^ to z = 1, which is much faster 
evolution than would be expected for a population usually 
found to have observational properties consistent with passively 
fading, old stellar populations. However, we cannot confirm 
this unusual result with our 0.75 < z < 1 .0 bin as we have just 
four galaxies in our sample and their mean surface brightness 
of 20.9 + 0.2 mag arcsec"^ suggests weaker evolution of only 
1.0 + 0.4 mag arcsec"^ to z = 1. A line of best fit anchored to 
the MGC value and passing between these two conflicting data 
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= 0.75 - 1.0 



2 = 1.0-1.25 
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Fig. 4. Galaxy morphology in our 'red, compact' class. Each row contains postage stamp images showing members of this class cut 
from the ACS/WFC /-band UDF image. These are clearly bulge-dominated systems, mainly giant ellipticals and SO/Sa galaxies. 
The rows are ordered in increasing redshift intervals from 0.5 < z < 0.75 to 1.25 < z < 1.5 to illustrate the changes in galaxy 
appearance due to bandpass shifting, decreasing angular sizes and morphological evolution. 



z = 0.75 - 1.0 



1.0-1.25 



z = 1.25 - 1.5 



Fig. 5. Galaxy morphology in our 'blue, compact' class. Each row contains postage stamp images showing members of this class cut 
from the ACS/WFC /-band UDF image. These are all small nucleated systems, featuring surrounding regions of irregular, diffuse 
material and only appear in our sample at high redshift. The rows are ordered in increasing redshift intervals from 0.75 < z < 1.0 to 
1 .25 < z < 1 .5 to illustrate the changes in galaxy appearance due to bandpass shifting, decreasing angular sizes and morphological 
evolution. 



points gives an evolution of 1 .65 + 0.22 mag arcsec ^ from z 
toz = 1. 



When moderated by inclusion of the 0.75 < z < 1.0 data as 
above our results are in good agreement with those of previous 
studies. For instance, Mcintosh et al. (2005) examined the 
luminosity-size relation of red sequence galaxies in GEMS. 
They found size-dependent evolution to z = 1 in V-band magni- 



tude of 1.6 + 0.1 mag for small systems (Re < 1 kpc), 1.3 + 0.1 
mag for medium systems (1 < < 2 kpc) and 0.7 + 0.1 for 
large systems (Rg > 2 kpc) (taking h - I for consistency with 
our work). In our 0.5 < z < 0.75 bin where we have enough 
galaxies to fit the luminosity-size relation slope we too find a 
suggestion of size-dependent evolution with a change in slope to 
-0.33 + 0.05 log(kpc)/mag from -0.281 ± 0.003 log(kpc)/mag 
at < z < 0.25, but this measurement is also consistent with 
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Fig. 6. Luminosity-size relations by galaxy type using our 'blue, diffuse' {top row) and 'red, compact' {bottom row) classes for 
the MGC 0.0 < z < 0.25 sample and the iUDF-BRIGHT volume-limited samples to z - 1.5. The two classes are seen to occupy 
distinctly different regions of this parameter space in the nearby universe (MGC sample). The 'red, compact' luminosity-size relation 
shows a higher mean surface brightness and a steeper slope than that for the 'blue, diffuse' type. These MGC relations are fit with 
contours and overlaid in grey on the high redshift iUDF-BRlGHT samples. Selection limits on object detection and reliability for 
the UDF survey are indicated with dark lines and grey shading. In the remaining whitespace areas galaxies are able to be detected 
with >90% completeness, and their sizes and magnitudes are able to be measured to within 10% of their true Kron values for more 
the 75% of cases. These boundaries were computed as described in Section 3, using the limit appropriate to « = 1 profile shapes 
for the diffuse types and that for n - 4 shapes for our compact types. The MGC selection boundary corresponding to an iso-volume 
contour at 100 Mpc^ is indicated with a thick black line. The bias-free regions inside the relevant selection limits of the two surveys 
are comparable for our volume-limited samples out to z ~ 1, at which point the observed distributions become overwhelmed by 
surface brightness-dependent biases. At z < 1 though, the UDF galaxies appear to show evidence for evolution towards higher mean 
surface brightnesses and we investigate this in Section |4!4l 



null evolution within the uncertainties. The conflict between the 
exceptionally strong surface brightness evolut ion detected in 
our 0.5 < z < 0.75 bin and the milder result of iMcIntosh et alj 
(l2005h (or indeed that of our 0.75 < z < 1 .0 bin) is presumably 
due to field to field differences (i.e., cosmic variance) across our 
small survey volumes. Catastrophic errors in the photometric 
redshifts computed for some of these systems may also be 
partially responsible. Finally, the 'red, compact' class suffers 
even more from the impact of selection effects at high redshift 
than does the 'blue, diffuse' class, so again we can gain little 
insight from our z > 1 samples. 



5. Discussion 

The bimodal nature of the bright galaxy population in vari- 
ous physical parameters is displayed locally in color-log(n) 
(Driver etal. 2006), log(«)-magnitude (Driver etal. 20()6; 
feall et all l2006i) and color-magnitude (iStrateva et alj \206)\i 
space , and at z ~ 1 in color-ma gnitud e space dBell et alj 
I2004t IMcIntosh etal] 120051: lBlantonll2006l) . This indicates the 
existence of two distinct galaxy sub-populations: (1) those 
of our 'red, compact' class that exhibit the red colors of old, 
passively-evolving stellar populations and the compact light 
profiles typical of spheroidal, velocity dispersion-supported 



structures, and (2) those of our 'blue, diffuse' class that exhibit 
blue colors indicative of recent star-formation and the diffuse 
light profiles typical of rotationally-supported systems. In 
gener al, as morphological comparisons indicate jDriver et al.l 
|2006| and here), these types may be understood to reflect the 
two sides of the Hubble sequence, namely ellipticals/lenticulars 
and spirals/irregulars. Driver et al. (2006) have argued that the 
bimodality specifically reflects the two component nature of 
galaxies with our two types corresponding to bulge-dominated 
and disc -dominated systems respectively. In any case, the 
existence of the galaxy bimodality suggests that there are two 
distinct pathways of formation, one leading to 'red, compact' 
types and one leading to 'blue, diffuse' types. Here we discuss 
these pathways and attempt to account for our results in light of 
contemporary galaxy formation theory and observation. 

Hierarchical clustering (HC) theory is the predominant 
galaxy formation scenario at present as it can link the evolution 
of the baryonic mat ter in galaxies to the evolution of their host 
dark matter haloes dWhite & Reeslll978|) inside the framew ork 
of the ACDM cosmological model dSpergel et al.1 l2007l) . In 
HC the two pathways of galaxy formation are explained in 
terms of the evolution of disc galaxies and ellipticals. Discs 
are thought to form from the collapse and subsequent cooling 
of primordial gas clouds inside rotating dark matter haloes 
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(e.g. iFall & Efstathloul [T980l) . spun up by mergers and the 
tidal torques of their neighbors in the early universe (^ Peeblesl 
[l969). Governato et al. (2007) find that in their ACDM-based, 
semi-analytical models a stellar disc forms in Milky Way-sized 
halos immediately after their last major merger, which is z > 2 
(as supported by the dark matter assembly history studies of 
iLi et al.l [20071) . At low redshifts disc evolution is characterized 
by accretion of small satellites, external gas accretion and 
secular processes such as bar formation. 

Elliptical galaxies are built up over cosmic time by m ergers 
between similar-sized discs jToomre & Toomr3 11972 ') and 
perhaps als o so-called 'dry-mergers ' amongst the e l liptica l 
population (Ivan Dokkumet al.1 1 19991: iNaab & Burker3 l2003h . 
Semi-analytical modeling confirms that the formation history 
of ellipticals in a A CDM universe follows that of dark matter 
halos dPe Lucia et al.l l2006t) with the most massive ellipticals 
and halos assembled at relatively low redshift. However, 
iDe Lucia et al.l (l2006l) find the bulk of stars in these systems 
were formed at high redshift in smaller galaxies, which were 
subsequently incorporated into the primary galaxy by merging. 
For instance, for ellipticals of mass >10" (Mb ~ -20 
mag) the median redshift when half of the stars are formed is 
~2.5, yet half of their stars are typically assembled in a single 
object only at redshift ~0.8. Hence, the galaxies in our sample 
(Mb < -18 mag) would be expected to show bimodality at z ~ 1 
in the HC scenario as both types are in place at this redshift, 
although the bright elliptical population should continue to be 
built up via merging to low redshift. This picture i s gene rally 
supporte d by observational studies such as Bell et al. (2004) and 
iBlantonI (2006) who find that the color-magnitude bimodality is 
in place by z ~ 1, but that the red sequence continues to evolve 
significantly in number density to the present day. In this study 
we have chosen not to examine evolution by type in the number 
density or relative fraction of galaxies in different classes as 
such measurements are especially prone to cosmic variance 
and selection biases. However, we can explore the pathways of 
galaxy evolution in terms of hierarchical clustering theory using 
the surface brightness and color evolution we have quantified 
for each type. 

For our 'blue, diffuse' (i.e., disc-like) systems we measured 
a reddening in (u - r),est color of 0.62 + 0. 10 mag and a decrease 
in surface brightness from z = 1 to z = of 1.57 + 0.22 
mag arcsec"^ in the B-band. The star-formation history of 
disc galaxies in the HC scenario is governed by the continual 
accretion of gas and dark matter into the parent ha lo, peaking at 
high re dshift but continuing on to the present dav. lWestera et al.l 
( I2OO2I) have used a HC motivated model to predict (B - y)rest 
color changes from ~0.30 mag at z = 1 to ~0.65 mag at z = for 
a disc at a 60 degree inclination angle. Using our SED template 
fits we convert our measured (u - r)iest colors to (B - y),est and 
find evolution from 0.36 + 0.10 at z = 1 to 0.59 + 0.22 mag 
at z = in good agreement with their results. We can also use 
this color evolution to compute a change in mass to light ratio 
(M/L), and thereby roughly disentangle luminosity changes 
from size chang es in the mean surface brightness evolution. 
Using iBell et all(l2003l) 's Table 7 we find our measured color 
change corresponds to a change in log(M/L) of 0.41 +0.11 
in the B-band, or a decrease in luminosity at fixed mass of 
1.03 + 0.28 mag. This implies that an increase in mean disc 
size of 22 + 14% from z = 1 to z = is required to account for 
the 1 .57 + 0.22 mag arcs e c'^ su rface brightness evolution we 
observe. ISomerville et al.l (l2008b have simulated the evolution 



of disc scale-lengths in a CDM universe. They find that although 
the host dark matter halos of discs continue to grow past z ~ 1 
the build up of halos from the inside-out means that the inner 
part of the halo where the baryons are concentrated changes 
little over that time. Their final prediction of 15-20% size 
evolution agrees well with our result. 



For our 'red, compact' type we measured a reddening in 
(u - r)rest color of 0.37 + 0.10 mag and a decrease in surface 
brightness from z = 1 to z = of 1 .65 ± 0.22 mag arcsec"^ in the 
B-band. Due to increases in computational power it has recently 
become possible to follow the evolution of disc and bulge 
sizes, as w ell as masses, in HC theory using semi-analytical 
techniques. lAlmeida et al.l (|2007) adopt an approach in which 
spheroid sizes are computed following a merger by applying 
conservation of energy and the virial theorem. These authors 
are able to reproduce many features of the local elliptical galaxy 
population, such as the Fundamental Plane, and offer specific 
predictions for evolution of the elliptical luminosity-size relation 
with redshift, which are illustrated in their Fig. 17. Essentially, 
they find luminous elliptical galaxies are lower in surface 
brightness by ~1 mag arcsec"^ in the B-band at the present 
day compared to z = 1, which is less than the 1.65 + 0.22 mag 
arcsec"^ we measure here. They note that their predicted surface 
brightness evolution is almost entirely due to a change in M/L 
of their passively-evolving stellar populations by a factor of 
~3 over this time. Using Table 7 of Bell et al.l (12003 ) again we 
can convert our measured color evolution to a M/L evolution, 
finding a change of Alog(M/L) - 0.26 + 0.11 (i.e., a factor of 
~1.8) between z = 1 and z = 0, or an decrease in luminosity at 
fixed mass of 0.65 + 0.28 mag. 



Recent evidence suggests that dry mergers may be an 
important component in gian t ellipt ic al galaxy evolution since 
z ~ 1 (Ivan Dokkumetal.1 fT999h . iKhochfar & STU (l2006h 
investigate the formation of elliptical galaxies in the HC 
scenario via semi-analytical modelling with a unique method 
of estimating post-rn e rger si zes, motivated by the finding of 
ISpringel & HernquistI (|2005|) that in elliptical galaxies the 
effective radius of stars formed during a merger is ~5.7 times 
smaller than the radius of stars that existed in the pre - merger 
systems (i.e., formed quiescently). IKhochfar & SiUd (|2006|) 
show that the scatter in their mo del merger-to-quie scent ratios 
behaves as a function of mass like lShen et al.l (l2003h 's measured 
scatter in local elliptical half light radii from the SDSS. They 
go on to predict mass-dependent evolution whereby the most 
massive galaxies evolve the most rapidly in size with redshift 
because they undergo more dry-mergers (and thus have higher 
quiescent star fractions) than do less massive ones. For galaxies 
in the mass range 5 x 10" Mq < M < 1 x 10'^ Mq (i.e., 
— 22.2 > Mb > -23.0 mag using M/L values based on our 
z = 1 colors) they find half light radii are ~55% larger now than 
at z = 1, while for galaxies in the mass range 1 x lO"* M© < M 
< 5 X 10" M0 (-18.0 > Mb > -22.2 mag) sizes are between 
~15 and ~25% larger. After subtracting the surface brightness 
evolution due to change in the M/L ratio we measured via the 
color change (i.e., 0.65 ± 0.28 mag) our results suggest elliptical 
size increases by 37 + 1 1 % since z - I - This is still faster than 
the predictions, but our results certainly favor a model with 
significant size and luminosity changes. 
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6. Conclusions 

We find that the color-log(«) plane can be a useful tool for 
subdividing the galaxy population in a physically-motivated 
manner beyond the nearby universe. We confirm that the 
locally-observed correlation between morphological type and 
color-log(n) plane position is in place by redshift unity. After 
separating both the UDF and MGC galaxy populations accord- 
ing to the saddle point in the (well-defined) local bimodality, 
we construct luminosity-size distributions by galaxy type in 
a series of volume-limited samples. By comparison of these 
distributions at high and low redshift we measure surface 
brightness decreases of 1.57 ± 0.22 and 1.65 + 0.22 mag 
arcsec"^ from z - 1 to z = for our 'blue, diffuse' and 'red, 
compact' classes respectively. These results are found to be in 
good agreement with predictions from recent galaxy formation 
simulations within the ACDM, hierarchical clustering paradigm. 
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